Measurements of ultracold neutron lifetimes in solid deuterium 
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We present the first measurements of the survival time of ultracold neutrons (UCNs) in solid 
deuterium (SD2). This critical parameter provides a fundamental limitation to the effectiveness of 
superthermal UCN sources that utilize solid ortho-deuterium as the source material. Superthermal 
UCN sources offer orders of magnitude improvement in the available densities of UCNs, and are of 
great importance to fundamental particle-physics experiments such as searches for a static electric 
dipole moment and lifetime measurements of the free neutron. These measurements are performed 
utilizing a SD2 source coupled to a spallation source of neutrons, providing a demonstration of UCN 
production in this geometry and permitting systematic studies of the influence of thermal up-scatter 
and contamination with para-deuterium on the UCN survival time. 

PACS numbers: 29.25.Dz, 28.20.Gd, 32.80.Pj 



Neutrons with kinetic energies less than 340 neV (cor- 
responding to a temperature T< 5mK) can be trapped in 
material bottles and are referred to as ultracold neutrons 
(UCNs) 0, |, §. UCN densities at reactor sources have 
gradually increased with reactor power and improved 
techniques for extracting the UCN flux. The highest 
bottled densities reported in the literature, 41/cm 3 , have 
been obtained at the Institut Laue-Langevin (ILL) reac- 
tor in Grenoble B. 

Measurements of the neutron electric dipole moment 
H [| and the neutron lifetime @, |8j D attest to the util- 
ity of bottled UCNs for fundamental experiments with 
neutrons. UCNs may prove useful in improved mea- 
surements of angular correlations in neutron beta-decay 
[flo| , |llf , although experiments of this kind utilizing UCNs 
have not yet been performed. All of these experimental 
programs have been limited by the available densities of 
UCNs. 

Superthermal UCN production, where the production 
rate of UCNs due to down-scattering in energy is larger 
than the combined up-scatter and nuclear-absorption 
rates in the material, was first proposed in 1975 by Golub 
and Pendlebury [Q in superfluid 4 Hc and experimentally 
investigated shortly thereafter |l3|, H. In this process 
phonon creation in the liquid is used to down-scatter cold 
neutrons to the UCN regime, while up-scattering is sup- 
pressed by maintaining the superfluid at sufficiently low 
temperature. Because He has no nuclear absorption, the 
only limitations to the density of UCNs accumulated are 
wall losses and neutron beta decay. The production of 



UCNs by this process has been observed and agrees with 
theoretical expectations @ H> @ 0. @- 

While superfluid 4 He is an excellent superthermal con- 
verter, a few other materials, such as solid deuterium 
(SD2), satisfy the criteria for superthermal production. 
The limiting UCN density, pucn, one can obtain using a 
SD2 source is given by the product of the rate of UCNs 
production in the solid, R, and the lifetime of UCNs in 
the solid, tsd- Pucn — Rt~sd- A storage bottle opened 
to such a source will come into density equilibrium with 
the density in the solid. This led to the proposal of a 
thin-film source where the inside of a neutron bottle is 
coated with a thin layer of SD2 and the bottle is embed- 
ded in a cold neutron flux j[9[ ^o). The volume comes 
into equilibrium with the UCN density with a time con- 
stant for the coupled system, r , given by, 
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where Vsd is the volume of SD2 in the source and V is 
the total volume of the storage bottle. A valuable sum- 
mary of SD2 thin- film sources is presented in The 
effects of gravity and of the potential of the solid as well 
as UCN losses other than absorption in the film have been 
neglected in this expression but do not fundamentally al- 
ter the above picture. Ultimately, the limit to the UCN 
density is established by the trade-off between the cold- 
neutron flux intensity and energy distribution (which de- 
termine the production rate) and the heat deposited by 
neutrons and gamma rays in the SD 2 and bottle walls, be- 
cause tsd is a strong function of temperature. However, 
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the predicted production rates in SD2 [M and lifetimes 
p2| have not yet been quantitatively verified. Efforts to 
utilize SD 2 sources at reactors identified possible gains 
but suffered from problems with cooling the solid at full 
reactor power [ p^[ . 

Pokotilovski pointed out the advantages of UCN pro- 
duction in SD2 at pulsed neutron sources J2j| and showed 
the UCN densities 4-5 orders of magnitude greater than 
existing reactor-based UCN sources might be possible. 
More recently, the use of spallation as a pulsed source has 
been suggested [^J], ^5| . In a spallation UCN source, cold 
neutrons are produced by moderating spallation neutrons 
produced in a heavy target by a medium-energy pulsed 
proton beam. These cold neutrons are used to drive a 
SD2 superthermal converter. In a spallation target, the 
amount of heating for each neutron is lower than in a re- 
actor, allowing higher neutron densities in the vicinity of 
a spallation target to be achieved. Even higher neutron 
densities can be ontained by pulsing the proton beam and 
valving off the UCN storage volume from the production 
volume when the beam is off and by using the time when 
the beam is off to remove heat from the deuterium. In 
this case, the maximum UCN density that is produced 
is limited only by the impulse heating of the SD2. Ex- 
periments with the stored UCNs can be performed while 
the beam is off, eliminating the backgrounds due to cap- 
ture gammas found near continuously operated reactor 
sources. 

In order to test the concept of a SD2-based spallation 
UCN source, we have built a test source that we have 
operated with single pulses of protons produced by the 
LANSCE 800 MeV proton accelerator at Los Alamos Na- 
tional Laboratory. In this letter we report the first mea- 
surements of t sd , the lifetime of UCNs in SD2 . Our mea- 
surements clearly demonstrate the critical influences of 
heating and para-deuterium contamination on the UCN 
lifetime, and provide a quantitative foundation for the 
development of SD2 superthermal sources. 

A schematic view of our apparatus is shown in Fig. [I]. 
Spallation neutrons were produced in a tungsten target 
with short (typically less than 160 ns long) pulses of 800 
MeV protons from the LANSCE accelerator. The fast 
neutron flux was contained and amplified using (n,2n) 
reactions in a layer of beryllium surrounding the spal- 
lation target. The spallation neutrons were moderated 
and cooled in a thin layer of polyethylene surrounding a 
58 Ni-coated stainless steel guide tube with an inner di- 
ameter of 7.8 cm. The polyethylene and the bottom of 
the guide were cooled with liquid helium, and a layer of 
deuterium was frozen on the inside of the guide. UCNs 
produced in the SD2 were confined by the guide tube and 
could be directed through a series of valves to the UCN 
detector. Neutrons were detected in a 5-cm-thick multi- 
wire chamber detector filled with a mixture of 3 He at 5 
mbar and CF4 at 1 bar. The low 3 He pressure and the 
large bend angle in the guide resulted in a high degree 
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FIG. 1: Schematic view of the apparatus used for this ex- 
periment. Neutrons were bottled in the region between valves 
A and B. Valve C was used to insert a thin 58 Ni-coated alu- 
minum foil in the guide in front of the detector. SD2 mea- 
surements were made by counting the number of UCNs that 
survive in contact with the deuterium as a function of time 
using valve A, with valve B and C opened. 



of selectivity for detecting UCNs in the apparatus. Data 
were acquired using a multi-scalar that was started by 
the proton beam passing through a toroidal pick-up coil 
and that scaled the count rate from the UCN detector. 

Several effects limit the lifetime of UCNs in SD2: up- 
scatter from phonons in the solid p0|, up-scatter from 
para-deuterium molecules in the solid 22| , absorption on 
deuterium, and absorption on hydrogen impurities in the 
solid. Model calculations exist for the contribution of all 
of these effects on the UCN lifetime in the solid. The 
total loss rate is a sum of contributions from each of the 
sources listed above: 
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+ l/Tpara + l/TDabs + l/THabs, (2) 

with the loss rate due to phonon up-scatter having differ- 
ent contributions from the ortho- and para-deuterium in 
the solid. Establishing the experimental basis to validate 
these models as both accurate and complete is essential 
for the design of a UCN source based on the superthermal 
production mechanism in solid ortho-deuterium. 

SD2 was frozen in the lower part of the cryostat us- 
ing a helium transfer refrigerator. The temperatures of 
the lower guide walls and of the liquid-helium cryostat 
were monitored with an array of silicon diodes mounted 
on the guide and the aluminum cryostat walls. The mea- 
sured temperatures tracked the vapor pressure curve of 
SD2 well at higher temperatures. The temperature of the 
solid at lower temperatures was obtained by averaging 
the temperatures of two diodes mounted on the outside 
of the guide wall. Later measurements made with diodes 
embedded in the solid indicate these measurements are 
accurate to I K. 

Both the hydrogen contamination and the para- 
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fraction in the SD2 were measured by means of rota- 
tional Raman spectroscopy on a gaseous sample taken 
by warming the deuterium after the measurement [ p2[ . 
These measurements yielded values for the HD concen- 
trations in the gas that varied from 0.2-0.3% with an 
uncertainty of about 0.1%. Other contamination was re- 
moved by a palladium membrane in the D2 gas system, 
before introducing D2 to the cryostat. The para- fraction 
was controlled by converting the D2 to a near thermal 
equilibrium otho/para ratio in an iron- hydroxide-filled 
cell p3] cooled to a temperature at or slightly below 
the triple point. In this way the para-fraction was re- 
duced from room-temperature equilibrium value of 33% 
to 2-4%. Intermediate values were obtained by mixing 
deuterium at room temperature equilibrium with con- 
verted deuterium before freezing. The precision of the 
para- fraction measurements was typically of order of 1%. 

The solid volume was measured by integrating the flow 
of gas while growing the solid. The volume was checked 
when the solid was warmed and the gas was returned 
to a buffer volume. The uncertainties in the pressure 
dependence of the calibration of the flow meter (20%), 
and uncertainties in the guide volume and temperature 
combined to lead to an uncertainty in the solid volume 
of about 20%. 

The sensitivity of the apparatus to UCNs was demon- 
strated by measuring neutron arrival times with and 
without a 58 Ni-coated 0.024 -cm-thick aluminum foil in 
place at location C (in Fig. |l|). These data are shown in 
Fig. |^. The number of counts, arriving in a time win- 
dow between 0.5 sec and 10 sec, with the foil closed was 
3% of the number with the foil opened. About half of 
these could be attributed to UCN leakage through the 
gap around the outside of the valve, and the rest were 
neutrons with normal velocities sufficient to penetrate 
the potential barrier provided by the 58 Ni. These data 
unambiguously demonstrate that the signal in the 3 He 
detector at the end of the guide was predominantly due 
to UCNs. 

If gravity, wall losses, the SD2 potential, and transport 
effects are neglected, and if the SD 2 is thin enough so 
that its volume is uniformly sampled by the neutrons, 
the lifetime of neutrons stored in a bottle in contact with 
SD 2 is given by equation ([!]). We have used this idea to 
measure the UCN lifetime in SD2. As depicted in Fig. 
[j], valves B and C were open for these measurements. 
UCNs were stored in the bottle between the end of the 
guide and valve A, which includes the SD 2 . The number 
of neutrons detected in a 10 sec wide time gate after the 
time, t, when valve A was opened was measured, typically 
for t =0.5, 1., 2., and 4 seconds. These data were fitted 
with the form c e~*/ T . The parameters c D and r were 
varied to produce the best fit for various conditions of 
the SD 2 . Background was subtracted by a fitting linear 
function to the data between 20 sec and 65 sec after the 
trigger and extrapolating under the data. The tsD were 
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FIG. 2: Background subtracted spectra with and without the 
58 Ni barrier in the beam at the location of valve C. The dashed 
curve is the result of a Monte Carlo calculation of UCN arrival 
times. 

extracted from the t using a lookup table of t$D vs t 
generated by fitting lifetime scans modeled with Monte 
Carlo for different tsD that were analyzed in the same 
fashion as the data. 

In the Monte Carlo transport, the full experimental 
geometry, gravitation, the SD 2 potential (108 neV), and 
wall collisions in the guide tube have been taken into ac- 
count. There are a number of parameters that must be 
determined from the experimental data in order to ex- 
tract tsD- the probability of UCN loss and nonspecular 
reflection for each collision with the guide wall, the SD2 
elastic scattering length, and the physical configuration 
of the SD2 frozen on the walls (SD2 can freeze as a flat 
"pancake" on the bottom of the guide, or coat the en- 
tire He-cooled surface at the bottom of the guide, in a 
"bucket" -shaped configuration) . 

Because we used very thin samples of SD2 for most of 
the lifetime measurements (0.4 cm), we were quite insen- 
sitive to the scattering length for UCNs in SD2 . The sat- 
uration of UCN production we observe for thicker sam- 
ples (up to ~6-cm of sample thickness) is consistent with 
the large scattering length calculated by Hill et al. p7[ . 
However, our extracted lifetime results change by less 
than 10% for scattering lengths as small as 0.5 cm. We 
therefore used the theoretical value of 8 cm for the elastic 
scattering length for all of the results presented here. 

The ratio of diffuse to specular reflections for wall col- 
lisions was adjusted to fit time-of-arrival spectra (see Fig. 

. The fraction of wall losses per wall collision and the 
shape of the SD2 ice were adjusted in a combined fit 
to reproduce the volume dependence for all of the low- 
temperature data. The extracted values for the guide 
parameters (0.025 for the ratio of diffuse to specular 
reflections) and wall losses (2.0 xlO -3 ) were consistent 
with those extracted from an independent set of mea- 
surements of guide transmissions and holding times at 
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FIG. 3: Left) Data points are measured SD2 lifetimes as 
a function of temperature, with the para-fraction fixed at 
2.5%. Solid lines show the predicted temperature dependence. 
Dashed line is the predicted effect of departure from the solid 
lifetime model due to up-scatter from the gas in the guide at 
higher temperatures. Right) SD2 lifetimes as a function of 
para-fraction for all of the data taken below 6K. The solid 
line is the model prediction of the para-fraction dependence 
at an average temperature of 5.6 K 



the ILL. The lifetimes due to absorption on deuterium 
and hydrogen were calculated using the known amount 
of hydrogen contamination, the tabulated thermal cross 
sections and the 1/v dependence of the cross sections. In 
addition the lifetimes in solid ortho- and para-dcutcrium 
due to temperature up-scattering were taken from the lit- 
erature 1 20, 2^]. The lifetime in solid para-deuterium due 
to molecular transitions was treated as a free parameter 
and found to be Tg°£ a = 1.2 ± 0.2 ms, consistent with a 
calculation that gives r|^ ra =1.5 ms [p2| . 

Results for UCN lifetimes tsd in SD2 as a function of 
the SD 2 temperature and para/ortho- fractions are shown 
in Fig. ||. The difference between the solid and dashed 
line demonstrates the the need to include the effect of 
deuterium vapor in the guide on the lifetime at higher 
temperatures. With this correction, the measured life- 
times agree well with theoretical predictions of the up- 
scatter rate. The main contributions to the UCN lifetime 
in SD2 have been measured and are quantitatively under- 
stood. These data demonstrate that UCN sources based 
on SD2 converters can provide the promised large gains 
over existing sources. 
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